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ABSTRACT 

SGR 1627—41 was discovered in 1998 after a single active episode which 
lasted ~ 6 weeks. We report here our monitoring results of the decay trend of 
the persistent X-ray luminosity of the source during the last 5 years. We find 
an initial temporal power law decay with index 0.47, reaching a plateau which 
is followed by a sharp (factor of ten) flux decline ~ 800 days after the source 
activation. The source spectrum is best described during the entire period by a 
single power law with high absorption (iVn = 9.0(7) x 10^^ cm~^); the spectral 
index, however, varies dramatically between 2.2 — 3.8 spanning the entire range 
for all known SGR sources. We discuss the cooling behavior of the neutron star 
assuming a deep crustal heating initiated by the burst activity of the source 
during 1998. 
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1. Introduction 

Soft gamma repeaters (SGRs) are a rare subclass of neutron stars characterized by their 
emission of randomly recurring outbursts of hard X- and soft 7-rays. There are currently four 
SGRs identified in our galaxy and one in the Large Magellanic Cloud (LMC), three of which 
have been found to pulse with periods ranging between 5 — 8 s (for a review see Kouveliotou 
2004). The very rapid secular increase of these periods, (spin down of P ~ 10""^° s/s), argues 
for angular momentum loss from a highly magnetized neutron star {B ~ 10^^ G). This idea 
was developed by Duncan and Thompson (1992, henceforth DT92) and, subsequently, by 
Paczynski (1992), following the detection of the most intense high energy transient observed 
to date, the giant flare of 1979 March 5 (from SGR 0526 — 66 in the LMC), which gave a 
lower limit of the source magnetic field of ~ 10^^ G. DT92 dubbed such sources "magnetars" . 

SGR 1627—41 was discovered with the Burst And Transient Source Experiment (BATSE) 
on the Compton Gamma- Ray Observatory (CGRO) in June 1998 (Kouveliotou et al. 1998), 
when it emitted over 100 bursts within an interval of 6 weeks (Woods et al. 1999); no further 
burst emission has been observed to date (July 2003). Roughly 20% of all events, represent- 
ing ~ 98% of the total burst-emitted energy (20 -150 keV), were bunched in an interval of 
three days, 1998 June 15-18. The X-ray counterpart to SGR 1627-41, SAX J1635.8-4736, 
was discovered in a BeppoSAX/Narrow Field Instrument (NFI) observation on 1998 August 
7, at a = 16^35"^49.8^ and 5 = -47°35'44" (J2000) with an error circle of radius V (95% 
confidence level; Woods et al. 1999). In Wachter et al. (2003) we report an improved (er- 
ror radius < 0.3") source location recently derived with the Chandra X-Ray Observatory, 
together with results of our near-infrared searches for a counterpart. A search for coherent 
pulsations in the BeppoSAX dataset, when the source intensity was higher, showed marginal 
evidence (~ 3a confidence level) near 6.4 s; to date the source spin period remains unknown. 

During their quiescent periods, SGRs have been identified as persistent X-ray sources 
with luminosities of ~ iQ^^sb ergs/cm^s. When active, their outbursts last anywhere from 
days to a year, with different burst frequency and intensity per source and per outburst. SGR 
bursts have typical durations of 0.1 s and their spectra are usually best-fitted to an optically 
thin thermal bremsstrahlung function with kT ~ 30 keV. The burst size distributions follow 
a power law of index —0.6, with peak luminosities of events ranging from the small, common 
bursts at <5xl0^^ erg/s to the rare, giant flares (only two have been observed) at ~ lO^^erg/s. 

The enhancing effect of the SGR bursting activity on the flux level of their persistent 
emission is well documented for SGR 1900 -|- 14 where the flux was shown to increase by a 
factor of ~700 after the giant flare of 1998 August 27. The flux from the source decayed 
within ~ 40 days according to a power-law, F oc (t—to)"^-'^ (Woods et al. 2001). For the next 
several months, the burst activity and persistent/pulsed flux level gradually declined. The 
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SGR was not observed to reach its quiescent flux level of 1 x 10~^^ erg/cm^s until 2000 April, 
approximately two years after the burst reactivation of the SGR, though extrapolation of the 
40-day decay after 1998 August 27 suggests that it may have done so between observations. 
The source has been intermittently active ever since, resulting in erratic flux decay behaviour. 

In contrast, SGR 1627—41 has remained dormant after a single active episode thus be- 
ing a better long term target for SGR cooling studies. We present here (Sections 2, 3, 4. 
and 5) the results of our monitoring of the flux decay of the quiescent X-ray counterpart of 
SGR 1627—41 obtained with imaging instruments (BeppoSAX, ASCA and Chandra) span- 
ning an interval of roughly five years. In Section 6 we discuss the source's atypical decay 
within the framework of neutron star crust cooling and its implications on the properties of 
SGR progenitors. 

2. BeppoSAX Observations 

Since 1998, we have observed the source four times with the BeppoSAX/NFl. Results 
from the first two observations were presented in Woods et al. (1999). Here, we present results 
from a refined analysis of the earlier data together with the two subsequent observations 
in August 1999 and September 2000. For all BeppoSAX data sets, we followed the data 
extraction procedure described in Woods et al. (1999). A log of the observations and spectral 
fit results is given in Table 2. Due to the signiflcant contribution to the background from the 
Galactic ridge, we arc using contemporaneous data from a concentric annulus around our 
source position for our background spectrum. An energy dependent multiplicative factor 
has been applied to the background spectrum to correct for reduced efficiency at off-axis 
angles in both the LEGS and MEGS. From these earlier BeppoSAX observations, where 
the higher source intensity allowed for better statistics, we have determined using XSPEG 
(vll.2; Arnaud 1996) that, in all cases, a single power-law (PL) model is the best fit spectral 
function for the data. 



3. ASCA Observations 

SGR 1627—41 was observed by the Advanced Satellite for Gosmology & Astrophysics 
[ASCA) on 1999 February 26-28. Onboard ASCA are four independent X-ray telescopes and 
four independent X-ray detectors; the latter are two Gas Imaging Spectrometers (GIS) and 
two Solid-State Imaging Spectrometers (SIS). For the present analysis, we have extracted all 
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available GIS and SIS data sets and applied the standard screening criteria^. We reduced 
the data separately for each GIS and for each SIS detector in 1-CCD mode, and extracted 
spectra using XSELECT (v2.1) from circular regions centered on the SGR with radii of 4' 
and 2.6' for GIS and SIS, respectively. Since the source was faint, we used a smaller region 
than the recommended ASCA source extraction region to reduce the background. Finally, 
the background was taken from neighboring regions on each detector. 

To derive an accurate measure of the source flux, we generated off axis GIS ancillary 
response files, dead-time corrected the GIS data, and included only the high and medium 
bit rate SIS data. Wc then generated response files (using FTOOLS v5.2 tasks as caarf and 
sisrmg) for each detector separately. We combined all the SIS data into a single data set, 
(similarly the GIS data), weighting in each case the response functions by their individual 
effective live times. We grouped the data into energy channels that contained at least 25 
events each. 

We performed independent fits to the SIS and GIS data (using XSPEC) of a PL model 
including the effects of Galactic absorption. We note that the final fiux we measure for the 
SGR (see also section 5 and Table 2) is significantly lower than the value reported earher 
{F — 5.1 X 10~^^ ergs/cm^ s) by Hurley et al. (2000). This discrepancy is due to the low 

source count rate in the GIS detectors, which results in a critical dependence of the measured 
flux on the selected background region. For example, choosing a background region several 
arcminutcs away from the SGR, where the detector count rate is much lower, wc find a flux 
comparable to the value reported by Hurley et al. (2000). In contrast, we have selected a 
background region close in angle to the SGR, which resides in an area of the sky that contains 
significant diffuse emission. This choice is justified by the fact that the fiux we measure with 
the GIS is in good agreement with the fiux measured by the SIS detectors (within 10%), 
which do not suffer from this issue due to their smaller field-of-view. 



4. Chandra observations 

We observed SGR 1627-41 with Chandra on 1999 September 20 (Obsl) and on 2003 
March 24 (Obs2). In both observations the SGR fell on AGIS-S3 (operating in Timed Ex- 
posure mode), a back-illuminated GGD with good spectral resolution. We used calibration 
products from Chandra GALDB (v2.21). Using the standard GIAO (v2.3) tools, we apphed 
the appropriate gains, CTI corrected the events, and filtered the data to include events 
with grades=0,2-4,6. We selected a circular extraction region of 3" radius centered on the 



^http: / /heasarc. gsfc.nasa.gov/ docs/asca/ screening, html 
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SGR for the source spectra. We collected the background spectrum from annular regions 
centered on the SGR with inner and outer radius of 3" and 30", respectively. The 2 — 10 
keV background-subtracted source count rates were 0.0049 cnt/s (Obsl) and 0.0058 cnt/s 
(Obs2) (the background contributed only ~ 3% to the total count rate for each observa- 
tion). We then generated the response files using standard ClAO tools and corrected the 
ancillary response files to account for the time dependent degredation of the AGIS quantum 
efficiency. We have investigated bremsstrahlung, blackbody, and PL models using XSPEG 
and find that the PL model gives the least G-statistic for both data sets (fitted individually 
or simultaneously). We have, therefore, adopted the absorbed PL model as the preferred 
continuum model for the Chandra observations. Our spectral fit results to the ungrouped 
data using the G statistic are given in Table 1. A simple ratio of the two data sets indicates 
that Obs2 has a softer spectrum, as also indicated by the spectra in Table 1. 

Table 1: Chandra/ AGIS Power-Law Model Fit Results 



Parameter 


OBSl 


0BS2 


0BS1+0BS2<^ 


A^h(10^^ atoms cm~^) 


9(2)" 


10(2) 


9.6(1.1) 


Index 


2.2(5) 


3.3(6) 


2.3(4), 3.1(4) 


Flux^ 


0.33 


0.44 


0.34, 0.39 


G-Stat (Goodness^) 


447.9 (86.6%) 


369.1 (92.6%) 


817.3 (95.2%) 



" Uncertainties are given at the 68% confidence level for 1 parameter ^ Unabsorbed 2-10 keV 

flux in units of 10^^^ ergs/cm^s 
Percentage of 1000 realizations of the model in which C-stat is less then that listed in table 
Obsl and Obs2 are simultaneously fitted with A^h forced to have the same value for both 

data sets and the remaining parameters are kept free 



5. Joint Spectral Analysis-Source energetics 

In Sections 2-4, each of the seven data sets was fit independently to a PL model atten- 
uated by interstellar absorption. Motivated by the consistency of the measured spectral pa- 
rameters from all seven observations, we performed a simultaneous fit to all data sets linking 
the Hydrogen column density, A^h- We have kept the value of A^h linked assuming that there 
is no obvious physical reason for the intrinsic source absorption to change while the source 
is in quiescence. For all observations, the fit was statistically acceptable {x^ /v = 690/655). 
We measure the (linked) effective hydrogen column to be Ah = (9.0 ± 0.7) x 10^^ cm~^. 
Table 2 fists all other fit parameters. 
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Table 2: Observation log and linked spectral fit parameters for SGR 1627—41. 



Mission 


Date 


Exposure 


Photon 


Flux" 




(MJD) 


(ksec) 


Index 


(10-12 ergs/i 


BeppoSAX 


51032.25 


44.9 


2.54(16) 


6.56 


BeppoSAX 


51072.52 


30.4 


2.77(21) 


4.94 


ASCA 


51235.20 


67.9 


3.24(24) 


2.76 


BeppoSAX 


51399.78 


80.4 


3.77(30) 


2.29 


BeppoSAX 


51793.59 


61.3 


2.99(29) 


1.86 


Chandra 


52182.50 


48.9 


2.17(30) 


0.267 


Chandra 


52505.00 


25.7 


2.95(36) 


0.266 



" Unabsorbed 2-10 keV flux 



Figure 1 displays the evolution of the source flux (upper panel) and spectral index 
(lower panel) since its activation. During the first ~800 days following the 1998 outburst, 
the fiux from SGR 1627—41 decayed monotonically as a power-law with exponent a = 0.47. 
The last BeppoSAX observation indicated that the source might have reached a plateau, 
approximately two years after activation. However, the two subsequent C/ianc/ra observations 
(at ~1200 and ~1500 days after the outburst) are significantly below (almost a factor of ten) 
the extrapolated decay trend. The photon index does show significant changes in the ~ 5 
years following the 1998 outburst; the probability that the observed variability is due to 
statistical fiuctuations is 7.1 x 10"^. We discuss these results further in the next section. 

While the spectra of SGRs 1900 + 14 and 1806 — 20 are relatively hard, with indices 
between 2.0 and 2.5, SGR 0526 — 66 is much softer with a spectral index of 3.5. Interestingly, 
the spectral index range of SGR 1627—41 (2.2 — 3.8) spans all other SGR sources observed 

at their quiescent states, although we do not sec the blackbody (bb) component observed in 
SGR 1900 + 14. This is not unexpected, given that the bb contribution in the unabsorbed 
flux of 1900+14 is of the order of 20%; since the column density along the line of sight for 
SGR 1627—41 is almost 10 times higher, extreme absorption between 0.1 - 2.0 keV would 
hinder the detection of any bb component that contributes less than at least 50% of the 
(unabsorbed) flux. 

Near the source activation the flux was 6.56 x lO^^^ ergs/cm^s, corresponding to a 
source luminosity of L = 9.5 x 10'^*^ crgs/s (assuming a source distance of 11 kpc; Corbel et 
al. 1999). The luminosity decayed slowly to 2.7 x 10^^ before it plummeted to its current 
value of 3.9 x 10^^ erg/s. In comparison, the persistent luminosity levels of SGRs 1900 + 14, 
1806 - 20 and 0526 - 66 are ~ 6 x 10^, 4.1x10=^^ and ~ 10^^ ergs/s, respectively. If 
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SGR 1627—41 decays further below its current luminosity level, it may well provide the first 
direct link between SGRs and Isolated Neutron Stars, whose luminosities it seemed to be 
approaching pretty rapidly through day 1200 (Figure 1). 

6. Interpretation and Discussion 

The study of the afterglow following an SGR episodic energy release is comphcated by 
the fact that both the magnetosphere and the surface are subject to change, and by the 
unpredictable nature of reheating by additional, later bursts. Conventional wisdom, though 
still somewhat tentative, would suggest that the surface radiation is reprocessed and non- 
thermalized by resonant scattering in the magnetosphere (Thompson, Lyutikov and Kulkarni 
2002) This suggests that sharp spectral changes would signify changes in the magnctosphcric 
configuration; conversely, an ordered, gradual decline in intensity with a constant spectrum 
indicates a decline in the thermal surface emission. To study the latter, a period of gradual 
decline is needed uninterupted by new events that reheat the surface. It was not obvious a 
priori that all these factors could be unraveled. 

The 40-day afterglow following the 1998 August 27 giant flare of SGR 1900 + 14 was 
extremely well fit by a crustal cooling model (Lyubarsky, Eichler and Thompson 2002, here- 
after LET), despite the existence of many small bursts that took place during that period. 
Ironically, the luminosity of this source was seen to increase and behave erratically over the 
following year, when there was relatively little SGR burst activity. In other words, SGR 
bursts do not necessarily heat the crust significantly and, moreover, there may be other 
kinds of sporadic heating that is not expressed in SGR bursts. 

Encouraged by the success of our outer crust coohng model for the 40-day afterglow 
of SGR 1900 + 14, we attempt here to understand the three year monotonic decline of 
SGR 1627—41 as cooling after a single deep crustal heating event coinciding with the burst 
activity of 1998. Details of our calculational methods can be found in LET. In particular, 
the leveling of the fiux during the third year followed by its sharp decline are curious features 
that beg for an explanation within this model. 

We present in Figure 2 our numerical calculations of the temperature evolution (cool- 
ing) of the neutron star crust with depth, assuming an initial energy injection to the 
crust of the order of 10"^^ ergs (estimates of the total energy released in bursts during the 
activation of SGR 1627—41 range between 4 x 10^^ — 2 x 10^^ ergs; here we assume that the 
conversion efficiency of the total energy released during the activation into soft gamma-rays 
is considerably less than 100%). The general shape of the initial temperature profile was 
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Fig. 1. — upper panel: Flux history of SGR 1627—41 derived with three different spacecraft 
over ~1500 days (2 — 10 keV). Errors in the flux are assumed to be 15%. The sohd hne is 
the numerical fit to the data as discussed in Section 6. lower panel: Evolution of the power 
law spectral index of the source during the same interval. 
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chosen under the assumption that the energy density of deposited heat varies with height 
more slowly than the specific heat: thus, in the outer crust, where the specific heat increases 
with depth, the initial temperature (i.e. immediately after the heating event) declines with 
depth. In the inner crust, on the other hand, the specific heat rises even more rapidly with 
the onset of neutron drip, but only if these neutrons are unpaired. (Paired neutrons do not 
represent new degrees of freedom because they are condensed into a superfluid.) Below Tc, 
the temperature below which neutrons pair, the specific heat of the crust material is small, 
and it rises dramatically near as free, unpaired neutrons appear. It is thus reasonable to 
choose the initial temperature to be near Tc in the lower crust. The theoretically predicted 
behavior of Tc as a function of density can vary from one model to the next, but the qualita- 
tive feature is that it increases with density and then decreases as the nature of the pairing 
changes, hence we put a bump in the initial temperature profile at 400 < z < 500 m. The 
exact shape of the bump is not very important as it gets washed out before it affects the 
surface temperature. 

We also assume that the core temperature is very low. This is justified if the mass of 
the neutron star exceeds 1.5Mq and the density, therefore, is high enough at the core to 
cool it via the direct URCA process. The assumption of a high mass neutron star makes 
the additional prediction of a thin crust, quite independently of the low core temperature, 
and this is reflected in the relatively rapid timescale (a few years) for cooling of the inner 
crust. If our interpretation of the transient cooling is correct, we thus determine the mass 
of SGR 1627—41 via two independent considerations to be above I.SMq. 

The time behavior of the surface luminosity is displayed in Figure 1 (upper panel), 
where we have plotted the data points together with the theoretical curve. It is seen that 
the plateau between days 400 and 800 is flt very accurately. The reason is that the sharp 
rise in the specific heat near neutron drip and above Tc makes this region in the depth 
temperature plane a very large heat reservoir, which keeps the temperature at neutron drip 
very stable. The duration of the plateau is determined primarily by the time required for 
the inner cooling wave to propagate outward to the neutron drip point. Most of the heat in 
the inner crust is conducted to the center of the star. The surface lightcurve is insensitive to 
many other details of the inner initial temperature profile, because by the time the surface 
feels the latter, the spatial variations within it are largely washed out. 

This model is unable, of course, to explain the March 2003 data point, which showed 
that the fiux did not decay further, and, moreover, showed a softening of the spectrum. One 
could argue that the apparent fiux "ledge" refiects a persistent, baseline luminosity that is 
observed only when the star is sufficiently cool. Other SGRs in fact have larger persistent 
luminosities, which are temporarily buried by enhanced, transient afterglows that typically 
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Fig. 2. — Evolution of temperature with depth, assuming a intial energy deposition of 10^^ 
ergs. The different temperature profiles correspond to years after the initial injection, as 
indicated in the plot insert. 
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follow major bursting episodes. In the other cases, (such as with SGR 1900 + 14), the 
persistent luminosity was attributed (LET) to a hot core (temperature ~ 7 x 10^ K). In the 
case of SGR 1627—41 however, the "observed" (based, for the time being, on only two data 
points) persistent emission cannot be attributed to a hot core according to our model: we 
calculated that any core hot enough to account for observable persistent emission at the level 
of the last two points in Figure 2, would keep the bottom of the crust warm and thus smear 
out the observed sharp ledge in the light curve. Further, the surface temperature would not 
have dropped fast enough to have accounted for the first Chandra data point. However, if 
the persistent emission is from a small hot spot at the surface, it does not affect the cooling 
of the deep crust and can be considered superimposed onto the theoretical cooling curve. If 
this low level component is ever found to be pulsed, it would support our conjecture that the 
baseline emission is heating at the surface rather than at finite crustal depth. More closely 
spaced and longer observations are certainly needed to establish the current level and nature 
of the source emission. 
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